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Abstract 

We study the behavior of the Binder cumulant related to long distance correlation 
functions of the discrete Gaussian model of disordered substrate crystalline surfaces. 
We exhibit numerical evidence that the non-Gaussian behavior in the low-T region 
persists on large length scales, in agreement with the broken phase being super-rough. 



1 Introduction 



The disordered discrete Gaussian model (ddg) is related to the 2D random phase Sine- 
Gordon model (rpsg). The two models belong to the same universality class, which 
includes different physical systems: for example it is supposed to describe crystalline sur- 
faces growing upon a disordered substrate [|l|, as well as 2D randomly pinned arrays of 
flux lines with the magnetic field parallel to the superconducting plane |@]. 

The standard approach used to investigate this model is renormalization group (RG) 
The application of the Mezard-Parisi variational approximation [^] (originally de- 
veloped for models with continuous replica symmetry breaking) leads to a one-step replica 
broken solution [§, and because of that it is problematic 0: still the main features found 
in the variational approach are very different from the ones found by using the standard 
RG. The most evident difference is maybe in the behavior one expects in the broken phase 
for the height-height correlation function. 



A considerable amount of work (see, for instance, 0-||T9[) has been devoted to the 
subject. Contrary to some former claims the situation has been shown to be in agreement 
with the existence of a broken super-rough phase, as implied by the RG approach, by the 
numerical simulations of p!9| . 

Recently very convincing further evidence for the system being super-rough down to 
T = has been established in references [^-|^, by exact computations of the ground 



state of finite volume samples of the ddg. 

In this letter we focus on the non-Gaussian behavior of the ddg in the broken phase. 
We show that it persists on large length scales. This is definitely not compatible with the 
predictions of the variational theoryQ On the contrary our data are well compatible with 
the behavior suggested by RG computations. 

Labeling with {di} the integer valued dynamical variables and by {rji} the quenched 
disorder, the Hamiltonian of the ddg model is 

^ ^ ^ (0, - <p^)\ ^.^d, + rj,, (1) 

<ij> 

where the sum runs over first neighboring sites of a bidimensional lattice. The model is 
related to the limit of coupling constant A ^ oo of the RPSG 

^^T.i'P^- ^3? - A E cos (27r(0, - r^O) , (2) 

<ij> i 

where now are the continuous dynamical variables. 
The partition function is defined as 

^.-Ee-^^ (3) 

where /3 = |t is the inverse temperature of the problem. 

A relevant observable quantity is the correlation function defined by 

^We remind again the reader that we are here in the framework of a one step rephca broken solution, 
and that this kind of criticisms Q does not apply to the variational theory when describing a continuum 
breaking ||]. 
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C(r,T) ^((0(ro)-0(ro + r))2). (4) 

In the high-T phase thermal fluctuations make the quenched disorder irrelevant. Both 
renormalization group and the variational theory predict a Gaussian behavior with a log- 
arithmic growth of the height-height correlation function 

CT>T.(r,T) --^ log r. (5) 

The critical temperature is expected in both approaches to be = ^. 

The renormalization group approach find that for T < one has a super-rough broken 
phase characterized by 

C^G^^(r,T)^alogr + MogV, 6 = (^Z£_I^' , (6) 

where a is a non- universal coefficient. 

On the other hand the Gaussian Ansatz of the variational approximation hints for no 
log^ r contribution, and the broken phase turns out to be described by a one-step replica 
symmetry broken solution with 



C^^^^(r,T)^^logr , (7) 



kn 

i.e. where the slope of the logarithmic term freezes at the critical point Tg. 

Let us recall in a few lines some of the main numerical results relevant for the prob 



lem. The authors of [|15[ studied RPSG without being able to detect any signature of the 



transition when measuring static quantities, probably ||T6| because of the small A value 



they used. Indeed for such small values of the coupling the difference from the pure case 
becomes sizable only on very large length scales. In numerical estimates for the cor- 
relation function C(r) of the DDG model were found to be compatible with the picture 
expected from the variational theory. Finally, evidence for the (logr)^ contribution to 
C(r) in the broken phase has been obtained in the case of the DDG model and in that 
of RPSG for different A values 



2 The Binder Parameter 

In this note we will mainly discuss about the distribution function of the height-height 
correlation functions at distance r. We define the probability distribution 



P[A(r), T] ^ {6 [A(r) - (0(ro) - 0(ro + r))]) , (8) 

where by (■ • •) we denote the thermal average (and here also an average over different 
values of fq), and by we denote the average over disorder. For sake of computational 
simplicity in the following we will only consider displacements r of the form (r, 0) or (0, r). 
The second moment of V is the usual height-height correlation function C(r), the focus of 
the investigation of J^, Here we will try to use the knowledge of the full probability 



distribution in order to gather more information about the system. 
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To characterize the probabihty distribution we measure 



D(r,T) = 3((A2(r)))'-(A4(r)) , (9) 
that together with C{r) allows us to define the Binder cumulant of V: 




mr)) 



(10) 



((A2(r)))V 

In the thermodynamic limit a value B = characterizes a Gaussian behavior. 

The work of [0 was based on the analysis of the Binder parameter for r = 1, that was 
providing evidence for a non Gaussian behavior in the low-T region. We will try here to 
answer some questions that are still open after |^ , looking at the long distance behavior of 
correlation functions: that will allow us to exhibit more evidence for a clear non-gaussian 
behavior. A purely Gaussian behavior could indeed be hidden for small r by short distance 
effects, and manifest itself only in the large r region. Analyzing B{r) in the large r region 
would make this effect clear. 

In order to understand what to expect for B we will use renormalization group (replica 
symmetric, to start with). We will mainly follow the approach described by Bernard in 
his Les Houches lecture notes |T^. In the field theoretical RG approach we start from the 



continuum version of (g) by writing 

5 = / ^ (f id,H^)f - A cos ($(x) - rf(x))) , (11) 

where the $(x) are the basic fields of the theory, and the d{x) are the quenched random 
field which make the system disordered. Universality is used to argue that these different 
systems exhibit the same critical behavior. Following Bernard renormalization calls for a 
generalization of this model: one introduces, in addition to the random phases, a random 
potential. If one would not do that at the start the random potential would in any case 
be generated by renormalization. So one writes 

, 2 

5 = / ^ i^,^x)f - A,ix)d,^x) - e(x)e^*(-) - r (a;)e-^*(-)) , (12) 
where the quenched fields ^ are distributed according to 

P[^] = e-^/'^^^* , (13) 

and the field can be written, by noticing that the rotational part decouples, as A^{x) = 
dfj_A{x), and it is distributed as 

P[A^] = e-^/'^^^^'^^' . (14) 

A U{1) symmetry guarantees to the model some remarkable properties. For example the 
(^-dependence of the correlation functions of the vertex operator e*°* can be factorized. 

Now one has to play the usual replica trick (at this stage with exact replica symmetry 
by definition). One writes an effective action for n replica's of the system, and takes the 
n — s> limit. The /3 functions of the model (for the running of k, g and a) turn out to be 
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(3^ = 2xa -2a'^ + ... , 

P. = 0, (15) 



where we have defined x = and n is not renormahzed. So in the low T phase, for 

K > Kc, it exists a non trivial infrared fixed point at cr*, with 

cr* = X + . . . for a; ^ 1 . (16) 

At cr* we have Pg = 9 ^^^^^ fiows (see Bernard [|13] for the characterization of such 

a run away fixed point) . 

We are interested in RG predictions for correlation functions. At the infrared fixed 
point one finds that 

G{r) = {e^iHr)~m))^ = r-2>e-4^(i°s")' , (17) 

where 7* is the anomalous dimension at the fixed point, 7^. = ^p* + 0{a^), and = 
1 + 0{x). We have already noticed that Pg^ = ^ (and it gives the large distance (logr)^ 
behavior of the correlation function). 

By expanding G in powers of a one finds that at all orders in perturbation theory (in 
the replica symmetric renormalization group approach) 



((<l>(r) - <l>(0))2) = ^ logr + ^(logr)^ , (18) 



and for small x 

^ x^ 



(($(r) - $(0))2) = - logr + — (logr)^ . (19) 



In the same way for the four point correlation function we find that 
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(($(r)-$(0))4) = 3M^logr + ^(logr)' -484 logr, (20) 

where we have defined d^, the unknown coefficient of the contribution to the anomalous 
dimension 7*, which we expect to depend from the temperature and that could even be 
zero. That means that in the RG approach we find that 

D^^^^(r,T)= 484 logr, (21) 
where D^^j,^ does not depend on /?*, and 



^RG . y^ ^ 244 logr ^ 244/^^ 

^<^=^' ' (l£.logr + ^(logr)2)2 log3 r • ^ ^ 
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So, if (i* 7^ 0, D{r) grows logarithmically with r. On the other hand according to the 
renormalization group picture linir^oo -B(r) = ~ (logr)~^. These are the theoretical 
predictions we use to interpret our numerical findings. 

As a last remark we want to notice that by studying a large version of the Random 
Phase Sine Gordon model, Bernard and Bauer ||12[ found (3g^, to be O (]^) • This means that 
there is no (logr)^ contribution to the height-height correlation function in the N oo 
limit in which the Gaussian Ansatz of the variational approach, corresponding to the 
leading order in a expansion, is expected to be exact. 



3 Numerical Results 



We have obtained our numerical data from simulations done on the Ape- 100 computer |^ . 
We have used square lattices of linear size L = 64 and L = 128, with periodic boundary 
conditions. We have fixed the surface tension k to 2. We have chosen the quenched random 
variables {rji} uniformly in the range (— |, +|]- 

We have simulated in parallel a total of 256 different realizations of the quenched 
substrate, and two uncoupled replicas for each sample. We have used a simple Monte 
Carlo local dynamics, by proposing to update in turn the {rfj} by an increment of ±1. We 
have used an annealing scheme, in which we have visited in turn decreasing values of the 
temperature T (T = 1.0, 0.95, 0.90, . . ., 0.40 for L = 64 and T = 0.90, 0.80, 0.70, 0.65, 
0.60, 0.45, 0.35 for L = 128). At each T values the thermalization sweeps were 0.5 million 
for L = 64 and 0.7 million for L = 128 (for further details see 0). 

We show in figure (1) B{r) as a function of T for different r > 4 values. We plot the 
data obtained on the smaller lattice, L = 64, since here we had a larger number of T values, 
but the behavior at L = 128 is very similar. When comparing L = 64 and L = 128 in the 
statistical precision of our runs there is no size dependence for r < 20. 

Deviations of B{r) from the Gaussian behavior at small distance become evident even 
for T < 0.85. The breakdown of the curves (for higher r values) in figure (1) is on 
the contrary compatible with the theoretical prediction Tc = ^. When looking at finite 
distances one finds a crossover (that does not correspond to a true critical behavior) for 



T > Tc [|T9|. Only measuring on large lattices real long distance properties one recovers the 
correct critical point, which turns out to coincide with good precision with the theoretical 
predictions. 

So, at lower T values the Binder parameter stays non-zero on larger length scales. What 
is even more important is the breaking of the slope of B versus T (see figure (1)). From 
our data we can say that for T < 0.65 the system is surely in its broken phase. 

In figure (2) we plot B{r) as a function of r for T = 0.45 and L = 128. Here the fast 
decay at short distance and the slow decay for large r is very clear. 

The data shown in figure (2) are qualitatively in very good agreement with the predic- 
tion of (p2D , with a non zero value of rf*. 

We show in figure (3) D{r) as a function of r for L = 64 at the highest temperature we 
have considered (T = 1.0) and at the lowest one (T = 0.4). That makes clear the difference 
between the high-T region, where D{r) is extremely small even at short distances and when 
increasing r becomes soon compatible with zero, and the broken phase where it is definitely 
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Figure 1: The Binder parameter as a function of T at different r > 4 values. L = 64. Lines 
are only meant to join neighboring points. Triangles (r = 4), squares (r = 5), pentagons 
(r = 6), hexagons (r = 7), etc. 




Figure 2: S as a function of r for T = 0.45 and L = 128. 



7 



0.015 




20 25 



Figure 3: D as a function of r at T = 0.40 (A) and T = 1.0 (□) for L = 64. 

non-zero and shows an evident increasing behavior. 

The fact that in the broken phase D increases with r is clear from figure (4), where 
D{r) is plotted as a function of r for T = 0.45, L = 128. That clearly shows that the 
non-Gaussian behavior of the model in the broken phase is not a short-distances effect. 

We have tried a quantitative analysis of the behavior of D{r), at a temperature well 
below the critical point. Following [0, [l^, we use in the fit the lattice Gaussian propagator 



1 1 - cos f 

- ^ „S „S 2 - COS (^) - M ^ N - ^ ^ . P3, 



ni=l n2=0 ^ - COS ^ j - COS ^ J 

which enables us to keep finite size effects under control. The data are very well fitted by 
the expected behavior D = ciP{r) +C2- Unfortunately errors grow quickly with r and data 
for r > 40 basically do not influence the fit. 

In figure (4) we show our best fit obtained by using data with 3 < r < 45 (disregarding 
ten more points at short distance does not change the results): 

lO^D = (114±3)P(r) + (36±3) , (24) 

where errors have been evaluated by using the jack-knife method. The value of the residual 
(per degree of freedom) is very good, close to 0.2 (but since the data points are very 
correlated the number does not have necessarily a deep meaning). The agreement with 
the renormalization group prediction (see equation (0)) is very good. Our best numerical 
estimate for is 7rci/24. 
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Figure 4: D as a function of r for T = 0.45, L = 128. The line is our best fit to the 
behavior D{r) = CiP(r) + C2 by using data with 3 < r < 45. 

4 Conclusions 

Our main conclusion is that the discrete Gaussian model for surfaces with a disordered 
substrate in the low-T region is non-Gaussian on large length scales. Such an evidence 
was needed to exclude the possibility of a short distance effect that could disappear in the 
asymptotic long distance regime. 

The picture which emerges from our analysis is therefore incompatible with the Gaus- 
sian variational Ansatz, calling for the broken phase being super-rough. Finally, our data 
give numerical evidence for rf* being non-zero, the behavior of D{r) at low temperatures re- 
sulting in good agreement with the logarithmic growth expected from the renormalization 
group approach. 
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